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Due to the attractive features that domain wall fermions possess with respect to chiral symmetry, we continue 
our investigation of the light quark masses with this discretization. Achieving reliable results, especially for 
(m„ + md)/2, requires strict control of systematic uncertainties. Our present results were obtained on a quenched 
P = 6.0 lattice with spatial volume « (1.5 fm)^. Consequently we remark on effects of finite volume as well 
as finite extent in the fictitious fifth dimension. We compute the renormalization factors nonperturbatively and 
compare to the one-loop perturbative result. 



1. INTRODUCTION 

Last year an exploratory calculation of the 
strange quark mass was completed using the 
domain wall fermion discretization within the 
quenched approximation |l[]. At the 15-20% 
level the results were encouraging: the pion mass 
squared extrapolated to zero at the input mass 
m — > 0, simulations at three different lattice spac- 
ings gave similar values for m^^(2 GeV) and the 
values were in agreement with other lattice re- 
sults. This work reports on the progress of the 
RIKEN/BNL/CU (RBC) Collaboration toward a 
precise calculation of the strange quark mass. 

2. SPECTRUM 

Using domain wall quarks we have computed 
light hadron masses on a 16^ x 32 quenched lattice 
with the following parameters: (3 — 6.0 (plaquette 
action), domain wall height M — 1.8, seven val- 
ues of valence quark masses (am's), and number 
of sites in the extra dimension Ng — 16. Results 
from 85 configurations for the pseudoscalar me- 
son, vector meson, and nuclcon masses are shown 
in Figures |] and || (see also Ref. ||]). 

With the light masses and improved statistics 
of this study compared to Ref. [Q , it has become 
clear that the linear extrapolation of (gMt^)'^ to 
am = gives a positive intercept, 0.018(2) (sta- 
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tistical error only), for these simulation parame- 
ters (see Fig. |^). Several effects could be respon- 
sible for this feature. First, the finite (and rela- 
tively small) volume of our lattice could increase 
the mass of the pion at a given am, and these 
effects would be especially visible at the lighter 
masses; as shown below the smallest am, 0.01, 
is roughly 1/4 the strange quark mass. Second, 
there is an intrinsic breaking of chiral symmetry 
due to the finite Ng'. the right- and left-handed 
surfaces states have some overlap within the fifth 
dimension and would result in a residual quark 
mass TTiros- Both finite V and Ng effects should 
be mostly am-independent. On the other hand 
studies by Columbia ||] at /3 = 5.7 suggest that 
in the large Ng limit one cannot account for the 
entire intercept by increasing the spatial vol- 
ume by a factor 2^. Quenching effects predicted 
by quenched chiral perturbation theory or due to 
artificial instantons could play a role. The issue 
of residual mass is addressed in Ref. ||] . 

We cannot presently distinguish am-dependent 
finite V and Ng effects from any nonleading or 
logarithmic behavior of the pseudoscalar meson 
mass. For this analysis we assume that the above 
effects are all am-independent, so then the chiral 
fimit is at am = —0.0038. We discuss systematic 
uncertainties regarding this assumption below. 

The am corresponding to the light quark mass 
is determined by extrapolating {Mt^/MpY to its 
physical ratio. There M^/Mp — 1.37(5) com- 
pared to the physical ratio 1.22. The lattice spac- 
ing is then set using the p mass, and the strange 
quark mass is set using either the physical K ot (j) 
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Figure 1. Pseudoscalar meson mass squared vs. 
am. 
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jeet to systematic error since disconnected graphs 
are being neglected, and so in the rest of this work 
we consider only ms{K). 

So far we have not completed an exhaustive 
study of the systematic errors at this lattice spac- 
ing. Therefore, for the present time we give rough 
estimates. Since Mf^/Mp is 10% higher than the 
physical value, our determination of the lattice 
spacing has an inherent uncertainty of at least 
10%. In this work we have assumed that finite 
V and Ng effects gave a positive m-independent 
contribution to the pion mass. If we instead as- 
sume such effects account for only half of the 
nonzero intercept ||] then the quark masses above 
would all decrease by 4 MeV. This is only a 5% 
change to the strange quark mass, but clearly 
the extrapolation to the average up and down 
quark mass is untrustworthy. A better under- 
standing of the systematic uncertainties is nec- 
essary before a reliable calculation of mi can be 
completed. The bare strange mass from this work 
is ms{f3 = 6.0) = 80 ± 6 (stat.) ±10 (sys.) MeV. 

3. NONPERTURBATIVE RENORMAL- 
IZATION 



mass. The bare quark masses (m^ = am-t-O.OOSS) 
are thus (with statistical errors only): 



bare 
mass 


1/a = 1.91(0.04) GeV 
lattice units MeV 


m 
ms{K) 
ms{4>) 


0.00166(0.00005) 3.17(0.11) 
0.042(0.003) 80(6) 
0.053(0.004) 101(7) 



The large difference in ms due to whether the 
K ov (j) meson is used to set the scale is a com- 
mon feature to quenched lattice simulations. This 
is succinctly expressed by the parameter J = 
MK-{dMv/dMlg) g. The value obtained using 
physical meson masses is 0.48±0.02 (the spread is 
due to the fact that (M^ - Mk'-)/{Mk- - Mp) ^ 
1); while our result is J = 0.34(7), comparable to 
calculations with Wilson fermions. In light of this 
difference, we do not average ms(K) and ms{(f>) 
but quote results for each input separately. The 
calculation of the (p mass in quenched QCD is sub- 



The mass renormalization constant has been 
computed perturbatively to one-loop for domain 
wall fermions 1^,^ . In this work we compute the 
renormalization constant nonperturbatively using 
the regularization independent momentum sub- 
traction (RI/MOM) scheme as advocated by the 
Rome-Southampton group ||]. 

To determine the renormalization constant for 
an operator O, we impose the following renormal- 
ization condition: 
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where Pq projects out the tree-level spin struc- 
ture of O, Ao is the amputated Greens function 
of O in momentum space, and Zq is the quark 
wavefunction renormalization. For fcrmion bilin- 
ear operators Ot{x) — q{x)Tq(x) only a single 
momentum space propagator is needed to com- 
pute the amputated Greens function 

Ar(p) = S-\p)(s{p)T{j,S\ph,))s-\p) (2) 



3 



where S{p) = exp(-ipj/)(q(?/)q(0)). In this 
work we compute momentum space propagators 
on 52 configurations 0. Eqns. (|l]) and (||) to- 
gether give Zr/Zq at any momentum p. In order 
for the RI/MOM method to work rehably the 
momentum should satisfy Aqcd <?C p <SC 1/a so 
that nonperturbative and discretization eff'ects, 
respectively, are small; however with presently 
accessible lattices this range is narrow. The ap- 
proach which we follow is to fit to the leading 
(ap)^ errors [B. First, we compute Zq through 
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which differs from Zq at 0{g'^) in Landau gauge 
H, then multiply Zr/Zq by Z'^ to give Z^ and 
divide by the two-loop running Cp M resulting in 
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which is scale independent through two-loop or- 
der. The result for the scalar density is shown in 
Fig. ||. The solid line is a linear fit to 0{o?p^) dis- 
cretization errors, and we take the extrapolated 
value at [apf = as our Z^'^^ . Zf^ at 2 GeV is 
obtained by multiplying by the two-loop running 
factor cfS(2 GeV). The result for = l/zf^ 
is 1.63(7)(9) (stat.)(sys.) which should be com- 
pared to 1.32, the one-loop perturbative Z,„ [|l],|j 
for these simulation parameters. 

Taking the lattice ms{K) from the previous 
section we find at /3 = 6.0 
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2 GeV) = 130 ± 11 ± 18 MeV 



(5) 



where the first error is the statistical uncertainty 
and the second is the systematic uncertainty. 

4. CONCLUSIONS 

We have completed a calculation of the light 
and strange quark masses at one lattice spacing 
and volume within the quenched lattice spacing. 
The increased precision has exposed finite Ng and 
V effects and these lead to errors of roughly 5% 
in the strange quark mass. We will soon have 
a calculation of on a coarser (/3 — 5.85) lat- 
tice, and of course further work at larger volumes. 
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Figure 3. Scalar density renormalization factor 
with perturbative scale dependence removed up 
to 0{g^). The solid line is a linear fit to a^p^ 
discretization effects. 



larger Ng, and smaller lattice spacing will reduce 
systematic uncertainties and yield a precise de- 
termination of the strange quark mass. 
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